Optical design of refiectionless complex media by finite embedded coordinate 
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Transformation optics offers an unconventional approach to the control of electromagnetic fields. 
A transformation optical structure is designed by first applying a form-invariant coordinate trans- 
form to Maxwell's equations, in which part of free space is distorted in some desired manner. The 
coordinate transformation is then applied to the permittivity and permeability tensors to yield the 
specification for a complex medium with desired functionality. The transformation optical structures 
proposed to date, such as electromagnetic "invisibility" cloaks and concentrators, are inherently re- 
fiectionless and leave the transmitted wave undisturbed. Here we expand the class of transformation 
optical structures by introducing finite, embedded coordinate transformations, which allow the elec- 
tromagnetic waves to be steered or focused. We apply the method to the design of several devices, 
including a parallel beam shifter and a beam splitter, both of which exhibit unusual electromagnetic 
behavior as confirmed by 2D full-wave simulations. The devices are designed to be refiectionless, in 
accordance with a straightforward topological criterion. 
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Metamaterials offer an enormous degree of freedom 
for manipulating electromagnetic fields, as independent 
and nearly arbitrary gradients can be introduced in the 
components of the effective permittivity and permeabil- 
ity tensors. In order to exploit such a high degree of 
freedom a viable method for the well-aimed design of 
complex materials is required. Pendry et al. reported 
a methodology based on continuous form-invariant co- 
ordinate transformations of Maxwell's equations which 
allows for the manipulation of electromagnetic fields in 
a previously unknown and unconventional fashion 
This method was successfully applied for the design 
and the experimental realization of an invisibility cloak 
Q. This work generated widespread interest specifi- 
cally in the prospects of electromagnetic cloaking, a topic 
that has dominated much of the subsequent discussion 
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transformation-optical designs reported in the literature 
so far have in common that the electromagnetic proper- 
ties of the incident waves are exclusively changed within 
the restricted region of the transformation-optical device. 
However, for continuous transformations the field manip- 
ulation cannot be transferred to another medium or free 
space and thus remains a local phenomenon. 

In this letter, we report a generalized approach to the 
method of form-invariant coordinate transformations of 
Maxwell's equations based on finite embedded coordinate 



transformations. The use of embedded transformations 
adds a significant amount of flexibility to the transfor- 
mation design of complex materials. For example, with 
finite-embedded transformations, it is possible to trans- 
fer field manipulations from the transformation-optical 
medium to a second medium, eliminating the require- 
ment that the transformation optical structure be invis- 
ible to an observer. The finite-embedded transformation 
thus significantly broadens the range of materials that 
can be designed to include device-type structures capable 
of focusing or steering electromagnetic waves. Like trans- 
formation optical devices, the finite-embedded transform 
structures can be refiectionless under conditions that we 
describe below. 

The mathematical formalism used for the calculation 
of the complex material properties is similar to the one 
reported in [|| [l6j]. Throughout this paper, we use the 
same terminology as introduced in these references. 

For a given coordinate transformation x a (x a ) — 
A" x a (j4" : Jacobi matrix, a = 1 ... 3), the electric per- 
mittivity t l 3 and the magnetic permeability /z 1 J of the 
resulting material can be calculated by 
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where det{(A\ ) denotes the determinant of the Jacobi 
matrix. For all the transformations carried out in this 
letter, the mathematical starting point is 3-dimensional 
euclidian space expressed in cartesian coordinates with 
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isotropic permittivities and permeabilities e y = EqS^ and 
fi lj — ^ S lj (S lJ : Kronecker delta) . 

A possible coordinate transformation for the design of 
a parallel beam shifter and a beam splitter consisting of 
a slab with thickness 2d and height 2h can be expressed 
by 



x'(x,y,z) 
y'{x,y,z) 



z'{x,y,z) 
with 
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where 2/ii is the maximum allowed width of the incoming 
beam, a determines the shift amount and I = 1 . . . n is 
the order of the nonlinearity of the transformation. 

The transformation equations are defined for (|x| < d), 
(\y\ < h) and (|z| < 00). For the case p = 1, equations 
©-© describe a parallel beam shifter whereas for p = 2 
the equations refer to a beam splitter. The Jacobi matrix 
of the transformation and its determinant are 
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FIG. 1: Linear spatial coordinate transformation of a parallel 
beam shifter. 



k[ := h(x',y') and 7' := 7(3/'). 

By equations ([I])-© it is straightforward to calculate 
the tensors of the transformed relative electric permit- 
tivity e r — e/eo and the relative magnetic permeability 
fi r = (j,/ no, which in the material representation are ob- 
tained as 
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where 
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is the metric tensor of the coordinate transformation. At 
this point it should be mentioned, that only the domain 
with 9(/ii— |/{|) = 1 has to be considered in the ma- 
terial implementation which simplifies the mathematical 
expressions. 

A linear coordinate transformation (1=1) for the de- 
sign of a parallel beam shifter is illustrated in Fig. [T] The 
space within the grey shadowed region I (<d(hi—\f[\) = 1) 
is tilted at an angle (f> — arctan(a) with respect to free 
space, which inherently results in a compression of space 
in region II and a dilution of space in region III. Consid- 
ering the boundary between region I and free space, the 
coordinate systems are seen to be discontinuous at the 
interface x — d. At this point, it is important to explain 
the difference between "embedded transformations" and 
" discontinuous transformations" . Interpreting the trans- 
formation as discontinuous, the boundary must be taken 
into account and the transformation of the y-coordinatc 
at the transition from region I and free space must read 

y'(x,y,z)=e(d-x)[y + aki{x,y)] +<d(x~d)y (18) 

so that a,2i\( x =d) 0(1 ${d — x). As in the ray approxi- 
mation the y-coordinatc lines in fig. [T] represent the di- 
rection of the power flow for a beam propagating from 
x = — 00 to 00. The inclusion of the delta distribution 
in the material parameters would result in the trajectory 
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FIG. 2: Nonlinear spatial coordinate transformation of second 
order for a beam divider. 



indicated by the black arrows. The incoming wave would 
be shifted in the y-direction and abruptly be forced back 
on its old path at the boundary at x = d, thereby ren- 
dering the entire beam shifting section invisible to an 
observer. However, as can be seen from equation ((4|, 
the boundary is not included in the calculation of the 
material properties for the beam shifter and beam split- 
ter. The coordinate transformation is carried out locally 
for the transformation-optical medium and then embed- 
ded into free space which results in the trajectory indi- 
cated by the green arrows. The beam is shifted in the 
y-direction and maintains its lateral shift after exiting the 
transformation-optical medium. This method is similar 
to the "embedded Green function" approach in the cal- 
culation of electron transport through interfaces [13] , so 
that we call it an " embedded coordinate transformation" . 

At this point, the question arises as to which condi- 
tions must hold for the embedded transformation in or- 
der to design a reflectionless optical device. We found as 
a necessary - and in our investigated cases also sufficient 
- topological condition that the metric in and normal to 
the interface between the transformation-optical medium 
and the non-transformed medium (in this case free space) 
must be continuous to the surrounding space. In the case 
of the beam shifter this means that the distances as mea- 
sured along the x-, y- and z-axis in the transformation- 
optical medium and free space must be equal along the 
boundary (x — d). As can be clearly seen from fig. [1] this 
condition is fulfilled by the embedded coordinate trans- 
formation of the beam shifter within the green shadowed 
region I. It should be mentioned, that the material prop- 
erties in the domains II and III can be arbitrarily chosen 
as no fields penetrate into these regions. 

As a second, more sophisticated example a beam split- 
ter is presented for the case of a nonlinear transforma- 
tion of second order. The material properties are de- 
scribed by equation (|16|) with (p = 2) and (1 = 2). This 
specific coordinate transformation is illustrated in fig. [2] 
The underlying metric describes the gradual opening of 



a wedge-shaped slit in the y-direction. The metric in the 
x-direction is not affected by the transformation. Similar 
to the parallel beam shifter, the beam splitter obeys the 
topological condition in order to operate without reflec- 
tion. 

To confirm our findings, 2D full- wave simulations were 
carried out to adequately predict the electromagnetic 
behavior of waves impinging on a beam shifter and a 
beam divider, respectively. The calculation domain was 
bounded by perfectly matched layers. The polarization 
of the plane waves was chosen to be perpendicular to the 
x-y plane. 

Fig. [3] depicts the spatial distribution of the real part of 
the transverse-electric phasor (color map) and the direc- 
tion of the power flow (grey lines) of propagating waves 
across a parallel beam shifter at perpendicular (fig. [3k-d) 
and oblique incidence (fig. [3b-f). The curvature of the 
incoming wave fronts was freely chosen to be plane (a-b) 
or convergent (c-f). As can be seen from fig. [3^-b, the 
beam shifter translates the incoming plane wave in the 
y-direction perpendicular to the propagation x-direction 
without altering the angle of the wave fronts. In contrast, 
the direction of the power flow changes by an angle </> = 
arctan(a) (a: shift parameter) with respect to the power 
flow of the incoming plane wave. After propagation 
through the complex transformation-optical medium, the 
wave fronts and the power flow possess the same direc- 
tion as the incoming beam, however the position of the 
wave is offset in the y-direction. The shift parameter a 
was arbitrarily chosen to be 1.8 (fig. [3^) and -1.8 (fig. 

Eb). 

A similar behavior can be observed for waves with 
wave fronts of arbitrary curvature, as for example for 
convergent waves (fig. [3]>f). In this case, the focus of 
the beam can be shifted within a plane parallel to the 
y-axis by variation of the shift parameter a. As is obvi- 
ous from fig.[3^-f, the same behavior applies for incoming 
waves at oblique incidence. The beam solely experiences 
a translation in the y-direction whereas the x-position of 
the focus remains unchanged. In all cases, the realized 
transformation-optical parallel beam shifter proves to op- 
erate without reflection confirming our metric criterion 
used for the design. The presented parallel beam shifter 
could play a crucial role in connection with tunable, re- 
configurable metamaterials as it would allow scanning of 
a beam focus along a flat surface without changing the 
plane of the focus and without introducing a beam tilt or 
aberrations. These properties become even more signif- 
icant for applications where short working distances are 
used between scanner and object. 

Similar simulations were carried out for the 
transformation-optical beam splitter. Fig. 0^ shows the 
electric field distribution and the power flow lines for 
waves at perpendicular incidence. The beam splitter 
shifts one half of the wave in the (+y)-direction and 
the second half in the (-y)-direction, thus splitting 
the wave at the mid-point. The split waves are not 
perfectly parallel at the exit plane of the device due 
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FIG. 3: Electric field distribution (color map) and power flow lines (grey) of a parallel beam shifter for diffracting plane waves 
with shift parameters (a) a = 1.8, (b) -1.8, for a convergent beam under perpendicular incidence with (c) a = 2, (d) = —2 and 
for oblique incidence with (e) a — 2 and (f) a = —1.2. 




Normalized Power Flow 1-100 Normalized Power Flow 1-100 



FIG. 4: Electric field distribution (color map), power flow lines (grey) (a+d) and power flow on a linear (b+e) and logarithmic 
scale (c+f) of a beam splitter for diffracting plane waves with shift parameters (a-c) \a\ = 15 for perpendicular incident waves 
and (d-f) |o| = 12 for oblique incidence. 



to diffraction of the incoming wave of finite width. As 
can be seen, there exists a small fraction of scattered 
fields within the split region which can be explained in 
terms of diffraction and scattering which is out of the 
scope of this letter. The beam splitter was found to 
operate without reflection in agreement with the metric 
criterion. 

Fig. [4Jd displays the normalized power flow inside and 
outside the device. In order to enhance the contrast in 
the visualization of the power distribution at the beam 
splitter output, the color scale is saturated inside the 
beam splitter medium. As obvious from the transforma- 
tion (fig. E]), the power density inside the transformation- 
optical medium is expected to be higher than outside the 
material, which is indicated by the density of the grid 
lines in fig. F2] and confirmed by the simulations. The 
power flow density abruptly decreases at the output facet 
of the beam splitter. For clarification, fig. F4J; shows the 
power flow on a logarithmic scale. By integration of the 
power density inside the gap region between the beams 
and the power density inside either the upper or lower 
arm of the split beams a power ratio of 10:1 was calcu- 
lated. The scattered waves in the gap carried about 4% 
of the total power. 

In fig. pjji-f the spatial distribution of the transverse 
electric field (d) and the power flow on a linear (e) and 



a logarithmic scale (f) are shown for an obliquely in- 
cident wave on the beam splitter. Again, the beam is 
clearly divided into two beams with a small fraction of 
diffracted and scattered fields inside the gap between the 
split beams. As for perpendicular incidence, no reflection 
was observed in the simulation at both the input and 
output facet of the beam splitter. Similar to the result 
for perpendicular incidence, the propagation directions 
of the outgoing waves are not parallel. In addition, the 
angles of the central wave vectors of the split beams with 
refer to the central wave vector of the incident beam are 
not equal. 

In conclusion, the concept of embedded coordinate 
transformations significantly expands the idea of the 
transformation-optical design of metamaterials which it- 
self was restricted to continuous coordinate transforma- 
tions so far. In contrast to continuous transformations, 
the expansion to embedded transformations allows for 
the first time to non-reversibly change the properties of 
electromagnetic waves in transformation media and to 
transmit the changed electromagnetic properties to free 
space or to a different medium in general. In order to 
design the medium as reflectionless, a new topological 
criterion for the embedded transformations was found, 
which imposes constraints to the metric of the spaces at 
the interface between the transformation-optical medium 
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and the surrounding space. This metric criterion was ap- 
plied in the conception of a parallel beam shifter and a 
beam splitter and confirmed in 2D full wave simulations. 
Both devices showed an extraordinary electromagnetic 
behavior which is not achievable by conventional materi- 
als. Both examples clearly state the significance of em- 



bedded coordinate transformations for the design of new 
electromagnetic elements with tunable, unconventional 
optical properties. 
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